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Abstract. It is presumed that the observed cosmic rays up to about 3 X lO'** eV are of Galactic origin, 
the particles being the ones which are found in the composition of the stellar winds of stars that explode 
as supernova into the interstellar medium (ISM) or into their winds. These pailicles are accelerated in the 
supernova shock. In order to obtain the observed cosmic ray spectrum it is necessary to take into account the 
diffusive losses in the Galaxy (which are making the energy spectrum more steeper). Another modification of 
the source spectrum is due to the fragmentation (spallation) of the cosmic ray particles, after their collision 
with the ISM atoms. In this paper we are proving that some particles are injected in the supernova shock one or 
two time ionized, and, also, that the contribution of massive stars (30Mo < M < 50Mo) accelerated particles to 
cosmic rays (where the winds are highly enriched in heavy elements) is 1:2 for elements with Z > 6. Another 
goal of this paper is to check if the particles are injected with the same velocity, energy or momentum. 



PACS numbers: 98.70.Sa 

1. Introduction 

Cosmic rays are particles, which have been accelerated in the Galaxy or into the extragalactic space. They 
come as electrons, protons, heavier nuclei, and their antiparticles. Into this radiation were detected almost 
all the elements of the periodic table as well as their isotopes. Only those particles are missing, that 
decay with a lifetime generally too short to be observed on Earth, like neutrons and many radioactive 
isotopes. The energy spectrum reaches from some MeV to 10^" eV/particle, the highest energy of any 
known radiation. From the lower end, up to the highest energies, the flux decreases by about 30 orders 
of magnitude. In consideration of such a large energy range, different detection techniques are necessary 
to explore the cosmic radiation. Methods used at medium and high energy (< 10'" eV) require detectors 
carried on balloons, rockets or satellites and, for studying the charged component of the cosmic radiation, 
this technique has been extended up to several lO'** eV/particle. At still higher energies the observations 
are done with detectors located deep underground and with air shower arrays, which cover the ultra high 
energy (UHE; > lO''* eV) and extremely high energy (EHE) regions (> lO'^ eV). Whereas the lower energy 
experiments at the top of the atmosphere measure the primary nuclei directly, the ground-based arrays only 
detect the air showers generated by the incoming particles. 

Up to now, the origin and acceleration processes of cosmic rays are still subjects of debate. As main 
sources of cosmic rays are assumed supernovae, quasars, active galactic nuclei (AGN) and Gamma Ray 
Bursts (GRB). Because they are deflected in the interstellar magnetic fields, the charged particles of the 
cosmic ray radiation appear nearly isotopic on earth, only the neutral constituents such as photons and 
neutrinos providing direct information on their acceleration site. 

The energy content of the cosmic ray radiation is about 1 eV/cm^ and, therefore, comparable with 
other forms of energy in our galaxy, like magnetic fields (~ 1 eV/cnr') or starlight (~ 0.44 eW/cm^). 

Recently were proposed three source sites for the observed cosmic rays ll2T]l : 

(i) Supernova explosions into the interstellar medium, or ISM-SNe. This component produces 
mostly hydrogen and the observed energetic electrons up to about 30 GeV, and dominates 
the all particle flux up to 10"* GeV. 
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(ii) Supernova explosion into predecessor stellar wind, or wind-SNe. This component produces 
the observed electrons above 30 GeV, helium and most heavier elements already from 
GeV particle energies. The component extends ultimately to several EeV. Since the winds 
of massive stars are enriched late in their life, this component shows a heavy element 
abundance which is strongly increased over the ISM. 

(iii) The radio galaxies produce a contribution which dominates beyond about 3 EeV, and consist 
mostly of hydrogen and helium, with only little addition of heavy elements below 50 EeV. 
At this energy the interaction with the microwave background cuts off the contribution from 
distant extragalactic sources: the Greisen-Zatsepin-Kuzmin (GZK) cutoff. There are a small 
number of events which appear to be beyond this energy, and whether they fit into such a 
picture is open at present. 

The theory was originally proposed in Biermann (1993, paper CR I) and in Rachen & Biermann (1993, 
paper UHE CR I). Various tests were performed in Biermann & Cassinelli 1993, paper CR II; Biermann & 
Strom 1993, paper CR III; Stanev et al. 1993, paper CR IV etc. 

2, Cosmic Ray Acceleration 

Cosmic rays with energies up to 100 TeV are thought to arise predominantly through shock acceleration 
by supernova remnants (SNR) in our Galaxy. A fraction of the accelerated cosmic rays should interact 
within the supernova remnant and produce gamma-rays. Recent observations above 100 MeV (the EGRET 
instrument on the Compton Gamma Ray Observatory) have found gamma ray signals associated with 
at least two supernova remnants - IC 443 and y Cygni. However, it is possible that the gamma ray 
emission from y Cygni to be associated with a pulsar within the remnant rather than the remnant itself. 
Further evidence for acceleration in SNR comes from the recent ASCA observation of non-thermal X-ray 
emission from SN 1006. Reynolds ifTTl and Mastichiadis lfT2l interpret the latter as synchrotron emission 
by electrons accelerated in the remnant up to energies as high as 100 TeV, although Donea and Biermann 
OJ suggest it may be bremsstrahlung from much lower energy electrons. 

Acceleration to somewhat higher energies than 100 TeV may be possible, but probably not high 
enough to explain the smooth extension of the spectrum to 1 EeV. Several explanations for the origin 
of the cosmic rays in this range have been suggested: reacceleration of the supernova component while 
still inside the remnant; by several supernovae exploding into a region evacuated by a pre-supernova star; 
by shock acceleration inside the strong winds of hot stars or of groups of hot stars. At 5 EeV the spectral 
slope changes, and there is evidence for a lightening in composition. Is likely that this marks a change from 
galactic cosmic rays to extragalactic cosmic rays. 

For stochastic particle acceleration by electric fields induced by motion of magnetic fields B, the rate 
of energy gain for relativistic particles with charge Ze can be written (in SI units): 



where ^ < 1 depends on the acceleration mechanism. 

The second order Fermi acceleration (Fermi's original theory) can be modified in the context of 
supernova (SN) shocks, or other strong astrophysical shocks, into the more efficient first order Fermi 
mechanism lfT6l . 

3. Propagation of Galactic Cosmic Rays 

3.1. The FIP Factor Correction 

The elemental and isotopic mass fractions of cosmic rays, when traced back to their sources, resemble 
general mass fractions (or with the general abundances - GA) of elements |19|, but display also some 
significant differences. The general mass fractions are based mainly on measurements of solar spectra and 
in carbonaceous chondrite meteorites. For the elements H, He and N the diff'erence between CR and GA 
are factors of 20-30 (underabundant in cosmic rays). The elements O, S, and Ar are underabundant in CR 




Cosmic Rays Origin and Propagation Model 



3 



by a factor of about 4 and the nuclide Ne by a factor of about 6. The elements C and Zn and the nuclide 
^^Ne are underabundant by a factor of about 2. 

The cosmic ray source abundances of elements and isotopes are deduced from the abundances of 
cosmic rays near Earth, applying propagation corrections for the secondary components produced by 
nuclear spallation in the interstellar gas. 

The cosmic ray source composition resembles that of normal stars like the Sun, modified by the 
photosphere, chromosphere to corona particle escape mechanism that result in a diminished element 
abundance for the elements whose first ionization potential (FIP) exceeds about 10 eV lfT9l . We note that 
FIP provides a partial organizing principle of the ratios of cosmic -ray source abundances relative to general 
abundances. The value of about 10 eV implies a first-stage injection temperature (at stellar photospheres) 
of ~ 10'* K (about 1 eV), with easier escape for the charged ones. When these particles become coronal 
particles, they are boosted up to energies near 1 keV. Some of the coronal particles are boosted up in 
energy, becoming flare particles, with energies of the order of 1 MeV. A further stage of injection (stellar 
flare particles near the shock waves of supernova remnants) at energies near 1 MeV is plausible, as relative 
ionization loss eff'ects on particle range and on the composition cancels due to the eff'ective charges of 
atoms near energies of 1 MeV/nucleon. 

It is to be noted that the elements (and isotopes of Ne) with FIP > 10 eV have a very wide spread, and 
that other organizing principles in addition to FIP are needed. Since cosmic rays' ^^Ne, C and O contain 
contributions from Wolf-Rayet stars, and Ne could be affected by suppression of light elements, it is useful 
to adopt the elements S and Ar for the FIP correction. A value of / = 0.27 for FIP > 10.4 eV is adopted, 
and I =\ for FIP < 8.4 eV. Zn, with FIP = 9.4 eV, is in the middle of the transition region. 

Other forms have been proposed for the FIP correction, a linear fit being still unsatisfactory (as seen 
in |[l9l) especially forN and ^^Ne. 

3.2. Contribution of Wolf-Rayet Stars to the Composition of Cosmic Rays 

A second difference between the cosmic ray source composition and the general abundance is the 
enhancement of ^^Ne, ^''Mg and ^^Mg. Prantzos et al. 1 15 1 showed that the abundance ratio C/O in cosmic 
rays (about twice the solar one) can be explained in terms of Wolf-Rayet star contribution to carbon. The 
winds of Wolf-Rayet stars are energetic, close to 0. 1 MeV/nucleon, relatively close to the energies of flare 
particles, and so, the acceleration of these wind particles is possible. The Wolf-Rayet stars go through 
two phases: WN, when the nitrogen produced in CNO hydrogen-burning cycle is abundant at the stellar 
surface, and WC, when the carbon, produced in helium burning, is abundant at the surface, and '*N burns 
into22Ne [,19J . 

3.3. Spallation Correction 

The spallation is one of the processes having the most important role in the modification of the abundance 
of cosmic rays in their transportation through ISM. 

The basic features of target fragmentation, sometimes called "spallation", are very well understood: 
heavy fragments arise from peripheral collisions of heavy ions or relativistic protons with the target nucleus. 
These so called "spectators" of the reaction are excited primary fragments which then decay into the final 
fragments by a sequence of evaporation steps. 

The spectrum of the residual nuclei seems to be determined to large extent, but not fully, by the 
evaporation process. 

The difference between the spallation in high-energy physics in accelerators and in astrophysics is 
that, in the first case, the targets are the nuclei and the projectiles are the high-energy protons, while in 
the second case the targets are the ISM protons (in first approximation) and the projectiles are the nuclei 
that form the cosmic radiation. So, everything depends just on which particle is putted into the reference 
frame. About 10% by number of the interstellar gas is helium. Hence, about 20% of cosmic ray generated 
nuclear spallation nuclear products are formed in nucleus-helium interactions. In the analysis of cosmic 
ray interactions with atmospheric nuclei (mainly nitrogen and oxygen) an accurate knowledge of nucleus- 
nucleus interactions becomes essential Il20l . 
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High-energy protons cause many different nuclear reactions and, in principle, all these nuclear 
processes have to be taken into account for a proper description of the total process. Possible outgoing 
particles are, for example, all light particles, gamma rays and (above incident energies of 150 MeV) pions. 
Furthermore, high-energy fission may occur and outgoing particles from the proton bombardment stage 
further reactions. Clearly, a wide spectrum of reaction products will be formed. Because of this diversity, 
the nuclear data needs for accelerator-based transmutation have been categorized in the following classes: 

(i) Proton-induced reaction data 

High-energy proton reactions represent the primary step of the total transmutation or spallation 
process, and the importance of nuclear data associated with these reactions is obvious. Also, secondary 
protons with lower energy than the bombarding energy are involved. Needed data include total, elastic, 
inelastic, spallation cross sections, fission cross sections and reaction cross sections for outgoing light 
particles (p, n, d, t, a, y, ji^, n^, ;r ) and fission products. In particular, total neutron and proton 
production cross sections are relevant. Energy-angle distributions are mainly required for outgoing 
high-energy neutrons and protons. Cross-sections for outgoing protons or neutrons constitute the most 
important data. For the method based on intense neutron fluxes, (p, xn) data have highest priority, but 
also the total proton, neutron and spallation yields are required. 

(ii) Neutron-induced reaction data 

After the interaction of a high-energy proton with a target nucleus, high-energy neutrons which have 
been knocked out by the intranuclear cascade will induce secondary nuclear reactions and/or fission 
products in the actinides and in the spallation target. The same data as needed for proton-induced 
reactions are required for these high-energy neutrons. The highly excited residual nucleus that is 
left after the intranuclear cascade also evaporates a large number of relatively low-energy neutrons. 
The low-energy reactions that are induced by these outgoing neutrons have been subject to extensive 
evaluations. The additional neutron-induced data that are required are therefore mainly those for 
energies between 20 and 1500 MeV. 

(iii) Other nuclear reaction data 

The spallation occurs when the supernova shock travels through the stellar wind and then hits the 
surrounding molecular shell of dense gas. Then, there are two ways to consider the spallation: first, in 
the approximation that the spallation is steady, and second, that the spallation takes place in an expanding 
medium with a time scale of expansion comparable to the time scale of the spallation. 

4. The Initial Mass Function (IMF) 

The IMF is assumed to be a time invariant mass spectrum with a power law of the form: 

n(m)dm oc m^^^^'^^'dm 

for which m is the mass of a main sequence (MS) star in units of Mq, mi < m < ot„ and n(m)dm - the 
number of stars in the mass interval m and m-\-dm; mi is a lower mass limit of stars, and m„ is an upper 
limit of stars. Normally, the IMF is derived from the observed present-day mass function (PDMF) in the 
solar neighborhood, which is assumed to be independent of time. The derivation of the IMF from PDMF 
is difficult, involving assumptions about the star formation rate during the lifetime of the Galaxy. For stars 
with lifetimes longer than the age of the Galaxy (m < IMq), IMF is derived by assuming an average star 
formation rate in the past, whereas for stars with lifetimes negligible relative to the age of the Galaxy 
(m > 2Mq), the IMF is derived by assuming a present-time formation rate and taking into account the 
stellar lifetimes (T(m)). Given the uncertainties in both theory and observation, the IMF variations can be 
parameterized, and the proposed IMF can be tested by means of a detailed chemical evolution model. This 
method has been adopted in many cases. For example, globular clusters are expected to have a IMF similar 
to the Miller-Scalo form, with a typical lower mass limit of 0.1 M©. The Miller-Scalo IMF can be well 
approximated by a half-Gaussian distribution in log m : 

0(log m) = Co exp [-Ci (log m - Ci)^] , 
where Q = 66.2, Ci = 1.15, C2 = -0.88 |10|. 
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5. Chemical Yields for MS and post-MS Stars 

To have acceleration of wind particles in supernova shocks, first it is necessary to have a star massive 
enough to make possible the supernova explosion and, also, in their pre-supernova stage, strong winds to 
enrich the circumstellar medium with heavy elements. 

In the Hertzsprung-Russell diagram, stars along the main sequence need to be considered in four 
separate zero age mass ranges: 

• Stars below about 8 Mq do not explode as supernovae; 

• Stars from about 8 to about 15 Mq explode as supernovae, but do not have a strong stellar wind, and 
so explode into the interstellar medium. In pre-supernova stage evolve as Red Supergiants (RSG); 

• Stars from about 15 to 30 have a substantial wind. This wind is enriched only in helium. The 
chemical composition of the wind at the time of explosion is approximately 50% in He and 50% in H. 
The mass in the shell of wind-swept material is moderate. In pre-supernova stage evolve as RSG; 

• Stars from about 30 Mq have strong winds. This wind is enriched in heavy elements and has little 
hydrogen left. The mass in the shell of wind-swept material is large. In pre-supernova stage evolve as 
Wolf-Rayet (WR) stars or as "luminous blue variables" (LBV) stars. 

The evolution of massive stars depends strongly on mass loss by stellar winds. The most illustrative 
example is given by the fact that massive stars with initial masses above 25 Mq are likely to be left with 
only about 5 Mq at the end of their Wolf-Rayet stage. Wolf-Rayet stars represent a late evolutionary stage 
of stars with masses of 25-50 which occurs as a result of the loss of all or at least the major part of their 
hydrogen-rich envelope. Although WR stars are among the most luminous stellar objects in the sky and 
contribute substantially to the chemical and dynamical evolution of galaxies, even their basic properties, 
e.g. their radii or their luminosities, are still not well known today. This arises mainly from the fact that 
WR stars are obscured by their own dense stellar wind, the origin of which is not yet understood and which 
renders the modeling of WR stars and their atmospheres difficult |6|. 

Through frequently modeled, the hydrogen and helium burning stages of massive stellar evolution 
are greatly affected by uncertain physics: convection, mass loss and the cross section for the ^^C{a, -y)'^0 
reaction. The extent of convective penetration into regions which are stable according to local criteria 
("convective overshooting") is especially uncertain. Observational arguments in favor of a moderate 
amount of overshooting is, for example, the main-sequence widening, but those arguments are indirect, and 
other processes such as rotationally induced mixing or improvements in radiative opacities might satisfy 
the same observational constraints. 

A second uncertainty concerns the efficiency of molecular weight gradients in preventing convective 
mixing ("semiconvection"). Two extreme assumptions are frequently seen in the literature: 

(i) that molecular weight barriers may be neglected altogether ("Schwarzschild criterion" for convection); 

(ii) that molecular weight barriers prevent almost any mixing ("Ledoux criterion" for convection). 

In between these extremes, semiconvection describes the rate of mixing in regions having 
superadiabatic temperature gradients, but stabilized by finite molecular weight gradients. 

The physics of mass loss is not sufficiently well understood for almost any evolutionary phase of 
massive stars to make a reliable quantitative prediction. Progress has been achieved recently in the theory 
of radiatively driven winds, which is applicable to the hot stages {Jejf > 15, 000 K) of massive hydrogen- 
rich stars, but this still involves considerable uncertainty and does not agree with observed mass-loss rates to 
better than a factor of 2. Much more uncertain are the mass-loss rates for red supergiants. Only approximate 
fits to the observations exist and a factor of 10 uncertainty may not be an overestimate. Finally, mass-loss 
rates for LB Vs and for the various types of WR stars are also uncertain. 

6. Mass Loss and its Effects on the Chemical Yields 

Mass loss by stellar winds is a dominant effect in the evolution of stars with an initial mass M > 20Mq. 
The standard picture for the pre-supernova structure of a massive star (at least near solar metallicity) is that 
of a red supergiant which contains a hydrogen-rich envelope of many solar masses at the time of explosion. 
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However, the inclusion of mass loss allows several alternative structures, all of which presumably have 
observational counterparts as supernovae. There may be red supergiants with a very small hydrogen-rich 
envelope, hydrogen-rich WR stars, helium stars (WR stars with a surface helium mass fraction Ys = 1 - Z, 
with Z as initial metallicity), or WR stars of WC or WO type, which contain large amounts of He, C and O 
(and no hydrogen) at their surface. 

The star may become a bare core, then enter the WR stage and follow all or part of the sequence WNL 
(L for late), WNE (E for early), WCL, WCE and WO. This sequence corresponds to a progression in the 
exposition of the nuclear products: CNO equilibrium with H present (WNL), CNO equilibrium without H 
(WNE), early visibility of the products of the 3a reaction (WCL) and the growing of the (C-HO)/He ratio 
in the WCE and WO stars. The comparison of the observed and predicted abundance ratios confirms the 
validity of the views on the CNO cycle and He-burning. 

The evolution of massive stars is highly dependent on the initial stellar metallicity Z through the effects 
of mass loss. The mass loss depends on metallicity because the radiation pressure on the envelope depends 
on the abundance of metallic ions. The mass loss M in O- and B-type stars varies like M ~ Z", with 
a value of a between 0.5 and 1 . When the star enters the WNE stage and after the disappearance of its 
hydrogen, a very important effect is that its mass loss rates essentially behave like M ~ M^^, probably 
with no dependence on the initial Z. Consequently, the M-rates are very large at the beginning of the 
WNE stage (up to lO^^Mo/year for M ^ SOMq) and they then decline very rapidly as the stellar mass is 
decreasing. This produces a remarkable mass convergence (for example, all stars at solar metallicity with 
initial M > 25Mq reach the pre-SN stage with a final mass of about SM©). The metallicity dependence of 
the mass loss rates has been included in various grids and models. As a result the relation between the final 
and initial masses is very different according to the initial Z. At solar Z one has the above-mentioned mass 
convergence, while at very low Z the final masses are very close to the initial ones. 

The initial metallicity of the star affects not only mass loss, but also the outcome of nucleosynthesis. 
For metallicities Z < Zo/20 mass loss has presumably a negligible effect on the yields of stars of all masses. 
During H-burning, the initial CNO transforms to ''*N, and part of the latter nucleus turns into ^^Ne during 
He-burning (through a captures and one /3 decay). '^C, ''^N and '^O all have equal numbers of neutrons 
and protons but not ^^Ne (10 protons and 12 neutrons). This surplus of neutrons (increasing with initial 
metallicity) affects the products of subsequent burning stages and, in particular, of explosive burning, fa- 
voring the production of odd nuclei ("odd-even" effect) lfT4ll . 

Mass Loss Rates 

• The mass loss rate employed during helium burning (up to carbon ignition) is M = -kM^-^ (with 
M in Mq and M in Mo/year) and k - 6 x lO^^Mo/year, so long as the carbon surface mass fraction 
does not exceed 0.02, and k - lO^^Mo/year afterward |8|. This is also the mass loss rate taken for 
hydrogen-free WR phase, but with k - 6x lO^^Mg/year |4|. 

• For WR stars with a nonvanishing surface hydrogen abundance (WNL stars) the mass loss might be 
assumed constant: 3 x lO^^Mo/year for M < 25 M^, 5 x lO^^Mo/yeai" for 25 Mq < M < 35 M^ and 
8 X lO^^Mo/year for M > 35Mq | 18 |. 

• For hydrogenless WR stars a mass loss rate as follows can be adopted: 



with a - 2.6, and A: = 5 x 10 ^Mo/year for hydrogenless WN stars (WNE stars) and ^ = 10 ^Mg/year 
for WC/WO stars Q. 

7. Effects of Rotation on the Surface Abundances and Chemical Yields 

For non-rotating stars, the surface enrichment in H and N occurs when the star reaches the RSG phase. 
There, CNO elements are dredged-up by deep convection. For rotating stars, N-excesses occur already 
during the MS phase, for solar metallicity Z = 0.02, the predicted excesses amount to factors 3 and 4 for 
initial v™, = 200 and 300 km/s respectively. At lower metallicity, the N-enrichment during the MS phase is 
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Table 1. Surface mass fractions of various isotopes in stellar evolution models in the initial mass range 
15Mo < Mini, S 50Mo at the pre-SN stage. The pre-SN configuration is also indicated, where RSG means red 
supergiant and WC stands for Wolf-Rayet star of the carbon sequence. The last column gives the initial mass 
fractions used in the stellar evolution calculations 1 9 |. 
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-4 


6.80 


10" 


-4 


6.80 


10" 


-4 


5.59 


10" 


-4 


3.63 


10" 


-4 


6.80 


10" 


-4 


2''Si 


3.44 


10" 


-5 


3.44 


10" 


-5 


3.44 


10" 


-5 


3.44 


10" 


-5 


1.06 


10" 


-4 


1.49 


10" 


-4 


3.44 


10" 


-5 


30Si 


2.28 


10" 


-5 


2.28 


10" 


-5 


2.28 


10" 


-5 


2.28 


10" 


-5 


7.86 


10" 


-5 


2.43 


10" 


-4 


2.28 


10" 


-5 


56Fe 


1.32 


10" 


-3 


1.32 


10" 


-3 


1.32 


10" 


-3 


1.32 


10" 


-3 


1.32 


10" 


-3 


1.32 


10" 


-3 


1.32 


10" 


-3 



smaller, due to the lower mass loss. However, there is a very large increase (up to a factor of ~ 10) for late 
B-type supergiants, because at low Z the star spends a lot of time in the blue phase and mixing processes 
have time to work. These predictions are in agreement with the observed excesses for galactic B- and 
A-type supergiants. Also, the very large excesses observed for A-type supergiants in the Small Magellanic 
Cloud are remarkably well accounted for 

The chemical yields are strongly modified by rotation. The larger He-cores obtained in rotating models 
at core collapse imply larger productions of He and other a-nuclei elements. This is the most important 
effect of rotation on the chemical yields. In addition, by enhancing the mass loss rates and making the 
formation of WR stars easier, rotation favors the enrichment of the ISM by stellar winds. 

The rotational diffusion during the H-burning phase enriches the outer layers in CNO processed 
elements. Some ''^N is extracted from the core and saved from further destruction. The same can be said 
for '^O and ^^Al, a radioisotope with a half-life of 0.72 Myr The mixing in the envelope of rotating stars 
also leads to a faster depletion of the temperature sensitive light isotopes, for instance lithium and boron 
||2]. The supergiants at low Z in the range of 9 to 20 Mq seem to be good candidates for the production of 
primary nitrogen IfTTI . 

8. Our Model 

After a discussion with P.L. Biermann we decided that, for computing the mass fractions of different 
elements in the winds of massive stars, it is useful to write: 

C Xi(M) <:>(M) M(Xi) dM 

Xi = -^^^ , (1) 

/^ ^(1)(M) mXi) dM 

where XjiM) is the mass fraction of element / (see Table [T]i in the stellar wind ejecta, for one computed 
stellar model (with given semiconvection and metallicity) and for one initial stellar mass. 
Here <1)(M) is the Miller-Scalo IMF: 
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<l)(logM) = Co exp [-Ci (log M - C2)'] , 

where Co = 66.2, Ci = 1.15, C2 - -0.88. M is the mass of MS star in units of Mq. 

In equation ([TJ, M,- and M/ are the initial and final masses between which the stars exploding as 
supernova enrich the ISM. For obtaining the contribution of all the masses from H-R diagram one must 
make a summation after all those intervals between 15 and 50 solar masses (so, for the stars that are having 
strong winds). 

The mass range is split in intervals because in a representation log Xj = /(log M): 

1) There is no unique function to describe the evolution (in logarithmic scale) of the mass fractions 
with the mass. In the first approximation, this evolution can be described by a step function. The first 
reason for this behavior is that stars from different mass ranges have different wind mass loss and so, 
different contributions to the enrichment of ISM. The second reason is that the mechanism through which 
the wind is driven changes at 25-30 M©; 

2) The composition in the wind between 15 and 25 Mq is dominated by He, and between 25 and 50 
Mo by C and O. 

The mass loss M, for the element / (in Mq) is: 

MiXi) = j Xi, ,urf(t, M) M(t, M) dt . 
Assuming that the mass loss is constant in time, the above equation will become: 

M{Xi) = Xi, ,„f{M) M(M) = Xi(M) M(M) ■ 

For practical reasons it is useful to put all these terms as function of In M instead of M. Then, the 
equation ([T]i can be written: 

Cm ^ Xf(M) <D(M) M M{M) dM 

~ i — M 5 

' Xm ^{M) M M{M) dM 

where we introduced the expression for M{Xi) and we assumed that X, ^ X, surf- 

We will have X,(ln M) - exp(T), where :F = a In M + ^ is the fitting function for hi X, = /(In M) for 
the domain in which the integral is defined. 



^(InM) = exp [in Co - Ci (0.434 x In M - Ca)^] . (3) 

M = exp(^), where = exp(lnA; + alnM), with k and a given into the section 6. @ can be taken as 
the logarithm of some numerical value of the assumed stellar evolution model, for the stellar mass range in 
which is done the integration and, most important, for the considered element. 

Introducing all this functions under the integral we can compute the mass fractions in the wind. 

In order to obtain the mass fractions in CR it is necessary to compute the ratios between the observed 
cosmic rays and the mass fractions from our proposed model for one example. Through this we will be 
able to evidentiate not only the general form, but also and some particular features of our formulas. 

Plotting the data of individual nuclei as provided by the direct experiments shows, in general, a fair 
correspondence in overlapping energy ranges. The spectra can be described by simple power laws in 
energy: 

<p = (poE^^ , 

where <p is the integral differential flux in particles/(m^ s sr TeV/nucleus), <po are the absolute flux 
normalizations and y the spectral index for various elements. The fitting was done by Wiebel-Sooth & 
Biermann. Thek results for elements between H and Ni can be seen in Table |2] 
In the particular case of, let say, carbon: 



(4) 
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Table 2. The spectral indices y and absolute flux nomialization (pQ for the various elements with nucleai' charge 
number Z. We can see that the spectral index for Hydrogen is close to 2.74, and for all the other elements 
(except Li, Be and B which are spallation products) it is close to 2.67 121]. 



Element 


Z 


00 [m^ s sr TeV/nucleus] ' 


r 


X 
df 


H 


1 


(10 57 + 30) 


- 10"^ 


2.76 ± 0.02 


0.70 


He 


2 


(6.73 + 0.20) ■ 


10-2 


2.63 ± 0.02 


2.10 


Li 


3 


(2.08 + 0.51) ■ 


10-3 


2.54 ± 0.09 


0.90 


Be 


4 


(4.74 + 0.48) ■ 


lo-'* 


2.75 ± 0.04 


0.37 


B 


5 


(-Q 95 + 79~) • 


lo-'* 


2.95 ± 0.05 


0.45 


C 


6 


(1 06 + 01) ■ 


10-2 


2.66 ± 0.02 


1.42 


N 


7 


(2 35 + 08) ■ 


10-3 


2.72 ± 0.05 


1.91 


O 


8 


(1.57 + 0.04)- 


10-2 


2.68 ± 0.03 


1.70 


F 


9 


(3.28 + 0.48) ■ 


10-4 


2.69 ± 0.08 


0.47 


Ne 


10 


(4.60 + 0.10)- 


10-3 


2.64 ± 0.03 


3.14 


Na 


11 


(7.54 + 0.33)- 


10-4 


2.66 + 0.04 


0.36 


Mg 


12 


(8.01 + 0.26) - 


10-3 


2.64 ± 0.04 


0.10 


Al 


13 


(1.15 + 0.15)- 


10-3 


2.66 ± 0.04 


1.24 


Si 


14 


(7.96 + 0.15)- 


10-3 


2.75 ± 0.04 


0.10 


Mn 


25 


(1.35 + 0.14)- 


10-3 


2.46 ± 0.22 


5.38 


Fe 


26 


(1.78 + 0.18)- 


10-2 


2.60 ± 0.09 


1.81 



where: 

ync = ync = JC ■ 

From (IDl: 

NciE) = N>ic(E) + Nuc(E) , 
which, in the general case, will behave like (j is for the isotopes): 

j 

with: 

".sum, I sum, i^^ ' ■ 

Then, taking the minimum energy at a given rigidity p* /Z and considering that is not existing one 
superior limit for energy at which the particle can be accelerated, the number of particles of one type will 
be: 

Now, because the mass fraction: 

Zk Nk nik ' 

with k = H,..,i,...,Fe , the observed cosmic ray mass fractions ratio will be (for carbon): 

Xc _ N^'' mc 

or, introducing A^^* and A^^* in the above formula: 
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XHe~ 1 c\'-^"' mne ' 

<p(\He r -^r- 

Pc P*H 

which can be simplified using the rigidity — = — - : 

X^^^ YHe-l IZcV'' IP^ilY"'''' mc_ 

Xhc (f>Q,He TC- 1 IZhJ \ £■() / niHe ' 

In equation (|6]l, - yjlniHeEi), with £0 = 100 MeV the lower energy cutofi^. 
The cosmic ray ratio between mass fractions will be: 

Xc Xnc + Xi3c 
Xhb Xhc 

with the general form: 

Xj _ YjjXij 

Xhc Xhc 

In our model, for obtaining the real CR mass fractions, is necessary to take into account also the FIP 
correction factors (for ionization losses): 

Xj _ TjjXj 

— ^element ^ 

Now, the mass fraction ratio of any element present in cosmic rays and He can be computed and 
compared with the model results for the same ratios. 

9. Results 

After several tests we saw that the stellar mass loss was just a second order correction, being possible to 
neglect it. Without this, the model mass fraction formula (O becomes: 

Xf(M) <:>(M) M dM 

i^'Xi(M) ^{M)MdM 

In our computations we had to consider also that, between two Table[T]mass values, the mass fractions 
in the wind are constant. 

First, we took an equal mass fraction contribution of RSG and WR stars to the production of the 
observed cosmic radiation. In this case, the obtained values are summarized in Table |3] where, in the 
second column is the nuclear charge number and in the third the mass number The column noted with FIP 
contains the first ionization potential for the considered element. In the ninth column we have the needed 
factor for equality between the model and cosmic ray abundances, each normalized to the corresponding 
mass fraction for helium, min and max are the minimal and, respectively, maximal factors with which we 
have to multiply the model mass fraction ratios in order to take account also of the spallation correction 
ratios and obtain the experimental value. 

\ model I \ observed i ^ \model 

min X I — '— < — '— < max X ' 



1 Ip*c 



c 



i-yc 



Xne/ \XHej \XHe 

In what follows we tried to prove that, after considering the spallation correction, the factor that 
differentiate the model mass fraction ratios from the observed ones (for even-Z elements) is related to the 
ionization losses (abundance decrease) for elements whose FIP exceeds about 10 eV, due to the particle 
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Table 3. Our model data considering that RSG and WR stars are having an equal contribution to the production 
of the observed cosmic radiation. 

(r, \model 
i FIP factor min max 

[eV] model/observed (for (for 



y observed ) y observed ) 



H 


1 


2 


1.051034 


0.0 




1.051034 




1.364058 


13.6 


1/1.42 


5.33 




He 


2 


4 


6.45098 


10-' 


1.31679 


10-1 


7.7052 ■ 10- 


1 


1 


24.4 








C 


6 


12 


4.22546 


10-3 


3.77525 


10-2 


4.19782- 10" 


-2 


0.05448 


11.27 


1/(1.42-43) 




2 


N 


7 


14 


7.45885 


10-3 


8.22535 


10-2 


8.97124- 10" 


-2 


0.11643 


14.52 


1/43 




2 


O 


8 


16 


1.6265 ■ 


10-2 


8.7725 ■ 


10-3 


2.50375 - 10" 


-2 


0.03249 


13.56 


1 




~ 1 


F 


9 


19 


5.99832 


10-^ 


5.47141 ■ 


10-12 


5.99837- 10" 


-7 


7.7848 - 10-^ 


17.4 


15-20 




2 


Ne 


10 


20 


3.29658 


10-3 


3.06398 


10-4 


3.60298 - 10" 


-3 


4.676 - 10-3 


21.48 


1 




~ 1 


Na 


11 


23 


1.11809 


10-* 


3.24637 


10-5 


1.44272- 10" 


-4 


1.8724- 10-* 


5.14 


1 


0.97 


1 


Mg 


12 


24 


1.04327 


10-3 


7.49583 


10-5 


1.11823- 10" 


-3 


1.4513- 10-3 


7.64 


1 


0.62 


1 


Al 


13 


27 


1.17598 


10-" 


1.41344 


10-5 


1.31733- 10" 


-4 


1.7096- 10-" 


5.98 


1 


0.55 


1 


Si 


14 


28 


1.31863 


10-3 


1.01951 


lo-* 


1.42058- 10" 


-3 


1.8436- 10-3 


8.15 


1 


1 


~ 1 


Fe 


26 


56 


2.5597 ■ 


10-3 


2.40742 


10-1 


2.80044- 10" 


-3 


3.6345 - 10-3 


7.90 


20 


1 


2 



escape mechanism from the photosphere and the chromosphere to corona, and, also, to the different 
radiative preacceleration (before the shock injection) of WR star wind particle species. Proving that the 
radiative acceleration plays an important role in the abundance of cosmic rays will be an independent study 
and, because of its complexity, will make the subject of another paper. 

Analyzing the data in the Table [3] we see that the column 4 mass fractions are too high. Second, 
the columns 10 and 1 1 factors are not even allowing the spallation correction for even-Z elements (see O, 
Ne,Si), and no FIP correction factor (see column 9) is evident for the elements with FIP> 10 eV. 

In this theoretical interpretation we can not count on a correct value for the observed cosmic ray 
hydrogen because its principal contribution to the interstellar medium enrichment and cosmic ray spectrum 
doesn't come from stars with \5Mq < M < 5QMq, but from other sources: low mass stars (ISM SNe and 
RSG SNe) and Active Galactic Nuclei (AGN) jets. 

So, it is necessary to consider that the WR and RSG stars are having a different contribution to an 
element abundance in CRs. The formula which describes this can be written: 



vmodel _ vmodel , / 1 ^ \ vmodel fQ\ 

where for the previous studied case we had a = 1/2. 

In the same time we need a formula for the radiative preacceleration correction in the stellar wind 
and which, after introduction of the spallative correction, to give an overabundance in the wind (for Mg, Si 
and Fe), relative to the observed CRs, situated between 4 and 6 (if we are taking 1 as reference level for 
CRJsurface abundances for C, O and Ne). 

For clarity, for a FIP correction "hidden" in the mass fraction ratio, observed/model it is the same as 
saying CR/surface. 

We took first the preacceleration correction as being: 



^.observed ^ /^xr„.-l ^ ^ ^^del 

In this way: 

^ \observed i ^ xy/jM-I 4 ^ -ytnodelJiSG . /1 \ -ymodel,WR 



(9) 



V I 1 V^, ] A. -\rmodel,RSG , /i \ ^modeLWR ' 



\ohserved 

- ' must be equal with the observed flux at 1 TeV: 
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Table 4. Mass fraction ratios for even-Z elements when Q" = | and the preacceleration correction factor is 



elem. 


Z 


/ - \ model 


/ \ observed 

(1:) 




factor 
observed/model 




X- 


H 


1 


1.4915 


0.105 ±0.075 




0.9936+ 1.1161 




C 


6 


0.03278 


0.035 ± 0.015 


0.876 ± 0.022 


1 


1 ±0.1 





8 


0.0293 


0.045 ± 0.025 


0.998 ± 0.0265 


1 


1.75 ±0.55 


Ne 


10 


4.8955 ■ 10-3 


(2.22 ± 1.57)- 10-2 


1.145 ±0.022 


1 


0.75 ± 0.02 


Mg 


12 


1.5336- 10-3 


> 5.8- 10-3 


1 


6.945 ± 2.936 




Si 


14 


1.9436-10-3 


(5 ± 2) - 10-3 


1 


4.6828 ± 2.02 




Fe 


26 


3.8032 - 10-3 


> 3.05 - 10-3 


1 


3.1156 ± 1.7623 





\ observed / , \ observed 

Xi \ / 00,1 



Xi,He I \ 00,1 

In the above equations y,-„j is the spectral index at injection into the SN shock. 

After computations it was clear that, for any value of the a parameter, the data obtained for min, max 
and model/observed were not in the limits of the physics hidden behind these factors {min > 0.5 and 1.6 < 
max < 2). 

None of our further tries, with similar forms of (|9]l, gave good results (even using jdif fusion = 7/3 or 
the Table|2]y values instead of y,„y = 8/3). 
The formula that worked was: 

-Tr \ observed i ry \ / v \ model 
Xj \ / Z, \ / A, 



Xfje) \ZHel \XHe 

where s is one parameter that was determined as being ~ 1 (see Table |4|i. The factor ( — '—\ seems to be 

related to a different preacceleration (and, consequently, a phase space separation) of the wind particles 
before their injection into the supernova shock of the WR predecessor. 

In dH), trying a - 2/3 (a contribution for RSG of 66% and for WR of 33%), we obtained for even-Z 
elements the Table|4]values. 

Using other a factors we obtained also good results for a = 3/4, a = 4/5 and c = 3/5 (see Tables |5] 
|6]and|7li. 



Table 5. Mass fraction ratios for even-Z elements when = j and the preacceleration con'ection factor is 
(^r) ■ '^'^^ columns that are missing from this table are similar with the ones from the Q' = | table. 



elem. 


Z 






S 


factor 
observed/model 


H 


1 


1.52547 






0.9715 + 1.0912 


C 


6 


0.02439 




1.041 ±0.022 


1 


O 


8 


0.02785 




1.022 + 0.0265 


1 


Ne 


10 


4.9329 ■ 10- 


-3 


1.142 + 0.022 


1 


Mg 


12 


1.5505- 10- 


-3 


1 


5.4399+ 1.4418 


Si 


14 


1.9632- 10- 


-3 


1 


3.6154 + 0.9561 


Fe 


26 


3.8316- 10- 


-3 


1 


2.2071 +0.8492 



We checked also for 1/3, 1/4 and 1 values, but the differences between the s values for the triplet C, O 
and Ne were to big to consider any of these as a valid candidate for the exponent of Z. 

Comparing the ratios observed/model from the Tables|4]|5]|6]and|2]it can be remarked that, for Mg, Si 
and Fe the values which enter in the domains of all three are for a = 2/3 {observed/model - 4.0088 + 4.878) 
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Table 6. Mass fraction ratios for even-Z elements when a = ; and the preacceleration correction factor is 
(^)*. The columns that are missing from this table are similar with the ones from the a = | table. 



elem. 


Z 


/ j^, \model 


s 


factor 
observed/model 


H 


1 


1.55016 




0.956 1.0738 


C 


6 


0.02015 


1.148 ±0.022 


1 


O 


8 


0.02722 


1.033 + 0.0265 


1 


Ne 


10 


4.9751 • 10-3 


1.138 ±0.022 


1 


Mg 


12 


1.5663- 10-3 


1 


5.2662 ± 1.3628 


Si 


14 


1.9824- 10-3 


1 


3.4966 ±0.9014 


Fe 


26 


3.864 - 10-3 


1 


2.1234 ±0.801 



Table 7. Mass fraction ratios for even-Z elements when ff = | and the preacceleration correction factor is 
i^-^)' ■ The columns that are missing from this table are similar with the ones from the a = | table. 



elem. 


Z 


/ ^, \model 


s 


factor 
observed/model 


H 


1 


1.43411 




1.0334-=- 1.1607 


C 


6 


0.04011 


0.764 ± 0.022 


1 


O 


8 


0.03017 


0.984 ± 0.0265 


1 


Ne 


10 


4.7768 - 10-3 


1.156 ±0.022 


1 


Mg 


12 


1.4917- 10-3 


1 


8.7264 ±4.7126 


Si 


14 


1.892-10-3 


1 


5.9638 ± 3.3042 


Fe 


26 


3.7116-10-3 


1 


4.2764 ± 2.937 
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Figure 4. CRs mass fractions to stellar surface mass fractions ratio, corrected for spallation and transport in 
the ISM, as function of FIP. Here a is 2/3. 



and for a - 3/5 (observed/model - 4.01385 ^ 7.2134). Taking C, O and Ne mass fractions as reference 
ratios for CR/surface and equal to 1, we can see that the behavior of this ratio versus FIP is having a very 
good resemblance with the ratio CR/GA described in "The FIP Factor Correction" section of this report. We 
should remember that it was theoretically predicted a same FIP factor for the elements Mg, Si and Fe. One 
confirmation of our a choices comes from the fact that it is known from other theoretical and observational 
sources that the hydrogen mass fraction is 3 - 4 times less than the total CRs hydrogen mass fraction (in 

\ model 

j is SO 

high. The FIP for hydrogen is 13.6 eV, making from it an element affected by ionization loss and placing 
it on the same underabundance (compared with the elements with FIP < 10 eV) level with C, O and Ne. 
But, in both a - 2/3 and a - 3/5 the observed/model factor for H is fitting the theory by being at least 3 or 
4 times smaller than the one for Mg, Si and Fe. 

Each one of our two a cases is having its good and bad sides. For a - Ijli we have for C, O and Ne 
the smallest difference in s, but the sum over all the mass fractions for this a is bigger than the one for a - 
3/5, this sum being closer to 1 for the last case (^ 1 .2379 for a - Ijli, and ^ 1 .1065 for a - 3/5). 

The final choice will be made by testing the ^^Ne/^"Ne and (^^Mg+^''Mg)/^'^Mg CRs ratios for both 
alphas. 

10. Test 

Neon is one of the most important elements in CRs because its two major isotopes, ^^Ne 
and ^°Ne, are formed by different processes of nucleosynthesis (^•'Na(p, a')^^Ne; '^©(ff, y)^"Ne; 
''*N(a, y)'^F(e^, v)'^0(a, y)^^Ne). One neon component measured in carbonaceous chondrites meteorites. 
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Figure 5. CRs mass fractions to stellar surface mass fractions ratio for Mg, Si and Fe, coirected for spallation 
and transport in the ISM, as function of a. 
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Figure 6. The observed ratio between CRs and He mass fraction as function of our model ratio between CRs 
and He mass fraction. 
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called neon- A, is believed to represent the primordial neon isotopic abundance at the formation of the solar 
system, and this value has been adopted as "standard" solar system composition: ^^Ne/^''Ne=0.122 + 0.006 

The abundance of ^^Ne/^*'Ne measured for the galactic CRs at the orbit of Earth is 0.54±0.07. Using 
standard solar modulation and others cosmic ray propagation models, the CR source ratio inferred is 
22Ne/20Ne=0.38±0.07 (20Ne/22Ne=2.6±0.5) |i3l. 

In our production and propagation of galactic CRs model, the ratio ^^Ne/^"Ne for q'=2/3 is equal to 
0.3777, in very good concordance with the above source ratio. For a-315 this ratio is 0.4772, larger than 
the CRs source ratio. We also found that (^^Mg+^^'MgV^^Mg^ 0.3076 (for q'=2/3) and 0.3208 (for a=3/5) 
which, compared with the solar system value of 0.27, suggests that the cosmic ray source and solar system 
material were synthesized under different conditions. The obtained values are also suggesting that a-2jl> 
is the right proportion which we were looking for 

11. Conclusions 

• The observed cosmic ray radiation at energy up to 3 x lO'^ eV can be generated by acceleration of 
stellar wind particles in supernova shock (see Fig. |6]l. 

• The contribution ratio of RSG stars and WR stars to the CRs mass fractions is 2:1. 

• Ionization loss is responsable for the underabundance in the observed CR elements with FIP> 10 eV. 

• We can consider Si as reference element for overabundance of elements with FIP < 10 eV relative 
to those affected by ionization loss (see Fig. |4]i. The FIP correction factor will be, in this case, 
4.4434 + 0.4346. In this way, we can see that the elements with FIP less than 10 eV are having an 
larger mass fraction than the elements with FIP greater than 10 eV (relative to Si), by a factor of ~ 4. 
This happens to be exactly Zj^,y^j.„o„, the initial degree of ionization squared. Therefore, cosmic ray 
particles of an element with an initial degree of ionization of two are more likely to be injected into 
the supernova shock. 

• The observed cosmic rays have an increased mass fraction by another factor of ■ ^'^i^ factor 

seems to be the result of a differentiated radiative preacceleration (phase space dispersion), before 
shock injection, of different elements found in the Wolf-Rayet stellar atmospheres. The reasons which 
are taking us to this conclusion will be detailed in another paper. 
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